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ABSTRACT 


Simple analytical methods for the determination of the angle of 
emergence from complex poles and zeros are developed. Several case 
studies of root loci emerging from complex zeros are analyzed to deter- 
mine criterion for the loci starting from complex pole to end on con- 
plex zero, as a pre-requisite study for "cancel-compensation" and self- 
adaptive compensators. 
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CHAPTER I 


INTRODUCTION 


Many closed loop transfer functions of control systems such as 
those representing structual resonances in missiles, fuel sloshing, 
shaft twisting, etc., have a pair of complex conjugate poles near the 
imaginary axis. If gain requirements are not too high, the effect of 
the roots near these poles is often neglected because other roots 
dominate. However, in the design of high gain systems, these roots 
near the imaginary axis may cross over into the right half s-plane and 
thus cause instability in the system. Figure 1 shows the effect of 
raising gain on the open-loop Bode diagram. 

Cancel compensation would be an effective solution to removing 
these complex poles near the jw-axis if it were not for (1) the diffi- 
culty in the design accuracy of these compensators, and (2) the poles 
themselves may not be fixed, and may move about in an area in the s- 
plane. 

The mere placing of complex zeros near the complex poles is no 
assurance that the root loci is connected between these poles and 
zeros. (See Fig. 2) When this loci does close on the adjacent com- 
plex zero, the root on this locus has such a small residue that it can 
be neglected when compared with the other roots. 

The question remains: What are the requirements for the location 
of complex zeros adjacent to complex poles for best cancel compensation? 


The search for an answer to this question initiated this study. Such 
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wErfect of a pair of Complex Poles near the Imaginary 
Axis on High Gain Systems 
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a location or area might make it possible to locate complex zeros for 
the design of self-adaptive compensators. 

First in this study, the angle of emergence from a complex pole 
and zero was studied to determine how it varied and how it was related 
to the transfer function. 

Next, 43 different case studies were programmed into the Control 
Data Corporation 1604 Computer to obtain 215 root locus plots. These 
plots and observations are contained in Chapter 3. 

The term "compensated system'' will be defined hereafter to mean a 
system with a pair of complex conjugate zeros relatively close to com- 
plex conjugate poles. The "uncompensated system'' will be defined as 
the system without the complex conjugate zeros, but will include the 


real poles associated with a "zero-compensator." 
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CME 22k Il 
ANGLES OF @USACeeCE 
1. Deterxination of the Ansle of Zuerngzence 
The angle of emergence of tne root locus from & complex 
bole of zero hes oeen defined in the literature as the angle 
of the tanyvent to the root-locus curve et the pole or zero as 
SOWwh. it ee Os AY test point 8, neer the com 2iex pole on 


’ 


the root locus vVisiles tne oGuaev.on 





ng a 


Waere W418 an OGd interme. or zero, £ 2e the tovel nuscer of 


cE 
zeros, M is the total rnuiitoer of voles, and -7 is the angele 

@l the Comolex D0le Galinitesigel ly Close to the west point. 
Fige 3 shows that 25 is measurea counter-clockwisce with the 
negative abscissa axis as reference, whereas the angle of 
emergence, Yo is measurec counterciocxwise with tne positive 
ebsciece exis as refevencée. 2% iS thuc seen what che enple 
of emer;ence ie equai to the an 


me irom the test polnt to 


the complex pole, using ooo vos resoonding references. Thus, 


y= [Py =? Yin aa (2-1) 


The last two terms in the above equation can be deter- 
mined witn the complex pole or zero (from whicn tne engle of 
emersence is to be detersined) es the origin of the co- 
ordinate systez, since the test point s, is intinitesizally 
close to the comolex sole (or zero). Thus, the orizin of the 
S-plane is translated to the desired comolex pole or zero by 


setting w= s +“-j4 where w is the new complex variable in 


al 








22.5 ‘ 
Angle of Emergence of complex pole, 


Angle of Emergence of complex zero, (4 
Angle from test point 8} to comolex pole, Le. 


i2 


the w-oplane. A test point W, near the new origin (pole or 
zero) which 16 on the root-locus must still satisfy the 

angle criterion (180 degrees or an odd multiole). Age wl 
aporoaches zero the sum of tne ansles on the root locus olus 
160 degrees would oe tne aenzsle froa the origin to the test 
S0int. When jwils identicelly zero, there is no areular con] 
tribution from the pole or zero &t the origin and thus the 
sum of the remaining angies is the aen@le of eueneence. For 
examuple, a characteristic equation is partitioned as follows: 


as om Ve) = \ 
' ~ ~ 
G BY —_ - a -* ce _ —_ 
( coe aa poeee eB = Be 





To find tne an-le of emer7ence at the complex onie located 


at -—2 +j5 in the s-plane, transfer the origin to this pole 






by Setting. w= ste -ye Me cs = we 2 +55 Then, 
(w—2+ 5542-5) (re 9 545-55) = -] 
Wrot gd) (wre 0419) (Hee 54 oe +a oe 5 == 
(weed) (y=1- 370) ee 
w—oti5) CweEsa5 7) (Cw) (wee ld) ~ w4 
(w+1) (wrl+gic) 
a (are) 





Define the angle of emergence from the comsolex oole, 


> = Ww ha P 0 (2-22) 
= sa. Mee ae oO. Ay. 2 See Le aed SS 
Ve = tan T tan I (180 + aan i) + van s 
+ 90) 
We= 0 + 64.5 ~ (160 + 112.6 + 32 + 90) 
Yer 30.9 


Now to find the angle of emergence from tne contslex zero 
located at -3 +j5 in the s-piane, transfer tne origin to 


this zero by Setting w= s +3 = 95 or S$ =w-= 3 + 35 
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(w— 3+ 7454+ 3-45) (w-34+55+3+j55) 
wW- 3735) (w-34+954+1 W-37)5+t2-95 


w)(w+3510) a 
W- 3435) (wt7+j5) (wel) (w-1+j10 
eng atl Cope eT Cred eek) (2-3) 


Define the angle of emergence from the complex zero, 






eA Seay (2-3a) 


(= 180 + tan) 4 + tan? ~ + tan + tan? a 
(= 180 a ee oe 8 
(Z= 161.5 


2e Complex Residues and the Angle of Emergence 

The closed loop transfer function can be expanded by 
partial fraction expansion and the residues evaluated by use 
of Heaviside's expansion theorem: (a) multiply function by 
(s - Dds and (b) let s = P+ giving the residue at Kx 
Using the previous example again, the residue angle at 
s = -2 +j5 is: 

(ee ee (s+#3-15)(s+3+)5) (2-4) 

re s(stl0)(s+2r+j5) ae 
It is seen that tnis residue evaluation is similar to eoua- 
tion (2-1) except for the -1 in the denominator which intro- 


-Auces 180 dezrees. Thus the angle of the residue at a com- 


plex pole differs from the angle of emergence by 180 degrees: 


Unn [kag + 280 (265) 


The characteristic equation of the preceding section 
may be partitioned as follows, where the zeros are now the 


Nooles" and the poles are now the "zeros": 
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The left side of this equation is seen to be the recisrocal 
of the closed-loop treneier 2-unction. Altvnover chilis te ean 
improoer fraction, the residues can still be evaluated by 


Heavieide's expansion theorem. Thus, the residue at 


9 = -3+j5 ls: 
[Ge Se 5 3555 coe 
ee. (3599) 43495419) 03+ 65 eee) (eee) | (2-7) 
7, ae Pad a=, 
K, = 341.5 


Cofoaring equations (2-7) end (2-3) 1% Ls geen that 
they differ by 180 dezrees. Tnus the anzle of emergence 
from a complex zero is 

= [K, + 160 (2-7a) 


wnere iS is determined as 1lllusetrateG above. 


3. Effect of Conzlex Zeros neer Cornlex Poles 

What effect does the relative position of comslex zeros 
adjacent to conplex poles have on the anzle of emergence 
from the conplex soles? To answer this question, firet con- 
sider the change of tne angle of emersence with tne zeros at 
constant raGius frou the »oles wniie varying the anruler 
DOSLUL In Of the zeroes wlth ressect to the comozex ocoles. 

It is assumed that the disolacement between complex 
poles and the distance between comozlex pole and its 
neighboring zero is of a ratio of approximately 10:1. As 


seen in Fig. 4, as the angle 8, the anzsular position of the 
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Effect of Zero Angular Position on Angle of Emergence 


336 


Complex zéro, increases, the engzgle of emergence increases 
in magnitude. The anrular contripution of { due to the 
lower half plane zero is negligible compared with © in the 
upper hnelf plane. Therefore e linear csiope of one exists 


between anguler position of zero and angle of emergence, or 


ou 
ae : (2-8) 


WW) 


Next consider Fig. 5 waich saows the effect of reaial 
position of compiex zeros (with consvaent ENneGler DOSition) 
on the angle of emergence. As tne zeros vary in rectel posi- 
tion, there is no anyuler chan.e cue to us9er half oiene zero. 
The chenge in lower nalf plane zero anzle is g. In addition 
to the assumption stated avnadve, the lin 2,25 is assumed 
perpendicular ta line Poh Yor siell cnanees in redial posi- 


tions ss  Tiven 


g= ae cos @ where dis the disolacement 


a-r sin 9 distance between complex poles, 
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Effect of Zero Radial Displacement on Angle of Emergence 
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Sine it was assumed that @/r = 10, ad” >») -2drsinO +r 


a2 _ cos 8 7 (220) 
OY a 


“re total deriviutive of the angle of emercernce is: 
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For small disol«acements, 
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adjacent coaplex poles and zeros, is allowed to oecone zero, 
or “perfect cancellation" with poles, tne anguler position 
of the zeros 18 undefinea. the constant Sis then evaluated 
as the sum of all remaining zeros in the system minus the 
sum of all remaininzs poles, weasured from the unper nalf 
plane "cancelled" coz2.lex sole. Thus in the example cited 
im Secon 1. is evaluated as 210.9. For 6 = 180, the 
angle of emergence Vo = 50.9 degrees. Tne seconc term on the 
right side of equation (2-10) is negligible conpared with 
the other terms unaer the afssumptions thet have been made. 
(See Fiz. 7) It is seen thet the angle © aiffers from the 
residue angle of the uncompensated system (no couplex zeros 
present) by 90 degrees. 
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The angle of emergence of the uncompensated system, 
Y% » is readily determined: 
= k 2-54 
Y Ko + 180 (2252) 


Y= 4 - 90 + 180 =4 + 90 (2-12) 


Equation (2-7) can now be written as: 


Y=0+ C988 rp. + Wo — 90 


IR 


Y“zar WY - 90 (2-13) 
Thus the angle of emergence changes directly with the 
zero oosition angle @. (See Fig. 8) 
The residue anzglesof the complex pole and zero are 


defined in equations (2-4 ) and (2-7) respectively and can 


be wWricten in.the following form:: 








[Ko = Joy Res (2014) 
By /? 5 /? 3 au DOLE Py 

/% = /P} (Po Paks (2-15) 
= at zero Zo 


The relationship between the pole residue angle, Noe: 
and the zero residue angle kK, Cen) now be written ac: 


hee * P20 + [25 + [Pg (2216) 


bey xX = (25 > and 4 = /> 1 The equation (2-16) becomes: 


hola bak a | (2-17) 


The angles < and 4 are. seen in Pigs6. 14 is observed from 
thas figure that “=°6 = 1760 “and | = 


Supeuitubine these relacionesapo imo cous.10m (2817) 


yields the following result: 


/X,, zs S/o ECO! aD (2-13) 
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Relationship between Residue Angle of Complex Pole 
and Couplex Zeroe ; 


Lig = Ly ta +7 


But from equations (2-5) and (2-7a), 
YZ - 180 


iL - 180 


bi 
i 


“ 


Therefore, 
Y- 180 = - “& + 180 + 20 - 180 or, 
G=-— 4% + 20 + 180 (2-19) 


Substituting equation (2-13) into equation (2-19) yields 
the relation between the angle of emergence from the com- 
plex zero and the angle of emergence from the complex pole 
when the complex zeros are not present. 

(= -0 -Y%+ 90 + 20 + 180 

Y= 0 -Y%- 90 (2-20) 
Using the relation from equation (2-5a), then 

Lx 9 -/k, + 90 (2-21) 
Subtracting’ «quation (2-13) from equation (2-20) gives 


the relation: 
4-Ke=-2Y= “2(0/K,, + 180) 


W~% = - fk \2n22y 


Thus the difference between the angle of emergence 


from a comolex pole and the angle of emergence of its 
nelghboring complex zero is & constant for a. given systen, 
and thus not a function of zero position 6. 

Fisures 7 and 8 are computer-solution graphs of the 
angle of emergence versus complex zero radial and angular 


positions respectively. 
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CHAPTER III 


ROOT LOCI LEAVING COMPLEX POLES 


Ge Data Designation 

Having analyzed the angles of emergence from complex poles and 
zeros, the next step in this study is to determine the effect of 
various system transfer functions on the root loci emerging from the 
complex poles, and to determine when this loci closes on its adjacent 
complex zero, i.e., the root loci starts at the complex pole and ends 
on the complex zero. As shown in Fig.2 of Chapter I, this loci may 
go to the zero at infinity instead of the complex zero. 

Over fifty different cases were tried and the data consisted of 
about 300 root locus plots using the computer program in Appendix I. 
Only 43 case studies (Fig. 9) are presented with 215 individual root 
locus plots (Fig. 10) for this analysis. 

Each case study or individual root locus plot is designated or 
classified using letters and numbers, such as "A-la'', where the first 
(capital) letter designates generally the number of real poles and 
zeros, and position of complex conjugate poles. The number indicates 
the variable position of the given number of real poles and zeros. 
The last small letter specifies the variable position of the complex 
zero with respect to the complex pole. 

Double, triple, etc., real poles are specified by placing an 'X" 


over an ''X"' or a group of "X's" near appropriate point on the plot. 
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Fig. 9d Case studies 
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The double real poles are arbitrarily used to represent the poles 
associated with a "'cancel-compensator" and is usually located, for 
this study, at -8 units except in a few cases they are placed at -4 
units. 

Case Ais a fifth order system. Case B is a fifth order system 
with one real zero. Case C is a sixth order system. Case D is a 
sixth order system with the complex conjugate poles moved out to -3 
+ j4. Case E is a sixth order system with one real zero. Case F is 
a seventh order system. Case G is a seventh order system with one 
real zero. Case H is also a seventh order system with one real zero. 
Case I is an eighth order system. Figure 9 presents these 43 case 
studies, and shows only the pole-zero configurations without the com- 
plex conjugate zeros. Figure 10 shows the root locus plots for these 
case studies when the complex zeros are placed in various positions. 

Table 1 is a summary of the graphical data from Fig. 10. In 
addition it contains the computed values of angles of emergence from 
pole and zero (for 9 = 0°) and the uncompensated residue angle using 
the formulas developed in Chapter II. A comparison with these calcu- 
lated values from the transfer function compares very good with com- 


puter root locus plots. 


30 








[ 


yee Ce = - 5 —e 





. 
—_ me ee ei ei ei 
‘ 





31 


Fig. 10a Root Locus Data 
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Fig. 10b Root Locus Data 


32 








| 

| 

| 

) 
QQ | oO 
\O \O 
I ) 
<i “1 

KY 

| 

| 

oxen ~ nee Se he ee 

5 | 

1 a 
i 

— 

l 

| 

\ 
.Q i oO 
iy “uy 
I, | ' 
<4 <s 

<K | 

| 

| 

ae 


{ 








0) 
\O 
I, 
<“h 
«> 
= _— te _ ~~ _—A _— _— ~— — —_— 
—_—_—_—_-—-S eo a. 
Oo 
LY 
t 
<4 
eo 


| 
| 
| 
| 
! 
| 
| 
| 
| 
] 
| 
| 
1 





Root Locus Data 


Figs 10c 


a2 





Fig. 10d Root Locus Data 
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Fig. 10e Root Locus Data 
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10£ Root Locus Data 


Fig. 


36 





Fig. 10g Root Locus Data 
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Fig. lOi Root Locus Data 
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Fig. 10k Root Locus Data 
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Fig. 101 Root Locus Data 


42 





‘ 
% 
—_—_—— = see ee i i ie 





Fig. 10m Root Locus Data 
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Fig. 100 Root Locus Data 
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Fig. lOp Root Locus Data 
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Fig. 10q Root Locus Data 





Root Locus Data 


Fig. 10r 


48 





Root Locus Data 


Fig. 10s 
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Fig. 10t 
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Fig. 10u Root Locus Data 


51 


Ze Observations from Data 

a. There is always a zero position angle 9 in the first 
quadrant with respect to the complex pole which will cause the root 
locus from the complex pole to end on the near complex zero. This is 
observed in all cases. 

b. If the root locus closes on the complex zero, then the 
locus from the pole moves in direction of least number of degrees 
(angle of emergence from pole compared with angle of emergence from 
zero). 

c. The locus closes on complex zero if there is a root locus 
segment on negative real axis to the right (closest to the origin) of 
the real part of the complex conjugate pole. (C-1, C-2, D-1, D-2, D-3) 

d. Just as the radial distance of zero had little effect on 
angle of emergence, it also has little effect on root locus plots ex- 
cept to make plots proportionately larger or smaller. (Compare Case 
A-1, A-5, and A-6). 

e. The locus from complex pole closes on complex zero if 
asymptote vector is 90 degrees and the asymptote centroid is far 
enough out on negative real axis. (i.e., Number of pole minus number 


of zeros is equal to two). See cases B-1l, B-2, B-3, and B-4. 
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TABLE 1 


SUMMARY OF DATA* 


Transrer Function: 







Other Data 


- - - - - = - -_ _ ~- - 
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A-3d | 270 NC 
Wee. SNe Dos ee Se 
| A-4a | 000 C 
A-4b 045 C 
| A-4c |090 ; C 
| A-4d =| 135 G 
|A-4e _—_-180 NC 
ia-Ge 2s §| te 
| Aste 1270 4% CC 
A=4h 4315 Cc 
tec? ee i * 
A-5a | 000 NC ' Same as A-1 except complex zeros are on .25 
A-5b =| (045 | C ; circle. 
1A-5c {090 } C 
i A-5Sd 1135 NC 
/A-Se 180 NC | 
WA-5e 48270 NC | 
wA-Se 2315 | NC 
= ae aa = Shai aa Mi es kata tena eae : 
| A-6a 000 | NC Same as A-l except complex zeros are on 1.0 
/A-6b = 045 Cc circle. | 
- A-6c 
A-6d 
A-6 





*"Desig" is the designation of each root locus plot, 9 (theta) is the 
angular position of the complex zero measured with respect to its 
adjacent complex pole, "'C'"' indicates that the root locus from the 
complex pole closes or ends on its adjacent complex zero, "NC" in- 
dicates that the root locus does not close on the complex zero. Values 
of /K_ were found by use of a spirule, and the values of GZ and 

were calculated by use of the formulas developed in Chapter II with 9 
equal to 000 degrees. 
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Table 1 (Cont'd) 





Kp = 170 +260 (= 2zO 





B-3a_ | 

B-3b 

B-3c 

B-3 | 

B-3 

B-3 | 

B-3g | 

ees 

B-4b- 

B-4 

B-4 

eg P (Geglawe lic 
B-5c [090 | | Ae te (Cate gt a+ (+44 
B-Sd | 225 | NC | 

B-Se | 270 | NC | 


CS FREEMAN cla 8 EO agree a: PA Chace A Dhan ton a EE lg A BT RR SERN A A Cece GABE TO RN NR IO aH tl ARCA te Ee A PD I AEN PR PE A IEE AO EPY 


ee ee % See 


ke =1a4 U£2alt YWY=326 
| | | 
: 
, 


“2 A+. URS 


ae ot poem See eS 
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Table 1 (Cont'd) 


Transfer Function: Other Data 


| 
a 


C 
C 
a Di Oe, el ee atl—- 4) 
Ca2e | 225 Cy J J 
c 
C 


tke #19 (2108 Ye=72 


(] 
i] 
i] 
ene! 
i] 
' 
* 


— 
fe 2) 
>) 





i) 
i) 
t 
' 
' 
i) 
p anes 
' 
t 
0 
' 


OURO te 0 enter so en - 





COG) Ca Cn Cy 


SIDR TOOTS ST APIECE A Bet 
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Table 1 (Cont'd) 





| Desig! 0 C,NC , Transfer poe Other Data 

eae | pe : | Ki (4+ a-- bata th) 

| D-le | 090 C 42 (242) (a+ g) 43474 Cetzgt 
D-ld 135 C 

eae |.“ fe ee 


leo iilas alas, 
So 
Qa Oo 
aS 
min © 
Aaa A aco 
Pa 
$ 
4 AW 
GS 
p(s 
& + 
f° 
rN co 
+ 
SER 
at fy 
yA 
rN + 
SiS 
at 
as 


| ie = See? “26s qZ= 7 





AQAaAaAaAaan 


rt Pa Aa 


Pr ag AES PR De Bt 3 A eee DI Epes Bel Be El Rin, Cte 


Sa ee OT Cee] eet ee Oe ee ere 





Table 1 (Cont'd) 





) Transter function: Other Data 


I< Cae) eee b(ata + ‘b) 
C 
C A(atd Ca+2) C24 lat a+lt4¥ 


ne | 7 Kp =H3s (2203 (4= 337 









Hesh. Gnoas) |) cc 
E-3c 090 NC 
E-3d !180 NC 
E-3e | 270 NC 
Ege hee a Koy 7/(-ata - L) ote ah 
E-4b ='045 G = ; a 
E-4c , 090 NC A (a+b) (a+ ey 4 +/+ 
E-4d 180 NC a | 
Ue s eat Wen sce @ 
| | ; 
se - ois ; 
‘ f 
| ! 
) | 
. 








t 


NC NiKe =9S YWstl8S Basse 


© 
t 
_ 
Qu. 
— 
[e.2) 
Oo 
Zz 
oe) 





Ey 


Desig 


G-2a 
G-2b 
G-2c 
G-2d 
G-2e 
G-3a 
G-3b 
G-3c 
G-3d 
G-3e 


Piper era) Lee es 


poeet ARR 


G-4a 
G-4b 
G-4c 
G-4d 
G-4e 
H-la 
H-1b 
H-le 
H-1ld 
H-le 
H-2a 
H-2b 
H-2c 
H-2d 


CO RO tn i i a I RE A I A A ta aN ss. DA) RE Aa Pe A 





H-3b 

H-3c | 
H-3d | 
H-3e | 


PIL bk LOR NO hcl nts ana Re Srl tle Pa te EAI, PE + 


SE DDE Ae NO SDA EOD 


,. -— @& ow 


EID o 


b 3 - - = 


H-2e [fh 270 


Table 1 (Cont'd) 


CNC Transfer Function; Other Data 
- | K (ats a+a~ ) a+& +- 4) 
c | 4 (a+3 <a+€)(42+8 “(ati-gé (2+! +44) 
ts [Kp = 50 (Lh =/40 (4 = 40 


@e-= se eF fF sw 28e#§ 829 & ££ = - se ese ses 282® esos FF 2 FF 2 2 2e F 83F 2 2 FF @& 


| aa 2 i ee (147) 
| 24 (+2 (Ae +3)" (a+ l-g4 atl tad 


/ Kp = 36 GZ 126 (= S# 


Zi Gi GiaGr GQ 


?) 





| allele a aad fatatyb) 
a (a+tit)(a-+3 AAS lat 3) (a4 galley 
NC [ip = (23 WZ =213 (Z= a 


«“_ es ses ese ese = ses ses ses 2s ss 2 @e@ 8s 2 8s 82« FF F ses 2wB ewe ese F F 2& 2F fF, @ 


oe) hea a = i Comenss a 
C (a+n(a en Bye! s)*(o+! —geXOst9) 


NC Le = =60 (2150 b= 30 


-— = = - = - -=— = - -— «= - - - -— = - - - -— = - - - —-— = -— - 


oe K (a+7) (a4+a- by (acted 

Cc 

C pater (246) ve (re 
NC 
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CHAPTER IV 


CONCLUSIONS 


Two analytical methods for obtaining the angle of emergence from 
a complex pole or zero have been presented. In the first of these 
methods, the origin was transferred to the complex pole or zero (w- 
plane) and then the angle of emergence was defined by equations (2- 
2a) and (2-3a). 

The second method was to determine the residue angle by means of 
Heaviside's expansion theorem and then simply add 180 degrees to the 
residue angle to obtain the angle of emergence. 


The following equations are exact, i.e., no approximations were 


Y = / Kpe i iz0 (2-5) 


(ZL = Kz + 18C (2-7a) 


used. 


The easiest method to obtain the angle of emergence is to find 
the residue angle with a spirule, then add 180 degrees. 


Other relations developed in Chapter 2 are: 


Gb =O +Ue-—70O (2-13) 


aa (L+26 +190 (2219) 

Y= o- Lh +70 (2-21) 

“== -2 ky (2-22) 
and since UY, =Lkp - 10 

(= 6 the + 70 (3-1) 


It must be remembered that these equations are only approximate 
with the linearizing assumptions and approximations which have been 
made. Figures 7 and 8 aré computer solutions of the angle of emergence 
vs. radial position of complex zero and the angle of emergence vs. zero 
positionangle respectively, and thus verifies the linear approximations 
for the assumed separation distance between complex conjugate poles. 

Thus knowing the angle of emergence from a complex pole YZ, the 
position angle ® of the neighboring complex zero, the angle of emer- 
gence from the complex zero is readily obtained by equation (2-19). 
Similarly if the angle of emergence from a complex pole of an un- 
compensated system is known, the effect of adding complex zeros on the 
angle of emergence is obtained by using equation (2-13). Equation (2- 
22) reveals that the difference between the angles of emergence from 
complex pole and zero is proportional to the "uncompensated" residue 
angle and is independent of zero position angle 9. 

The angle of emergence from complex pole and zero help to pre- 
dict the manner in which the loci leave the pole or zero. A vector 
is drawn from the pole in direction of the angle of emergence and a 
second vector is drawn from the zero in the direction of its angle 
of emergence. Then the tails of these vectors are extended until 
they intersect at a point and thus defines an angle less than or equal 
to 180 degrees. It was observed that whenever the root locus from a 
pole closes on its adjacent zero, the loci moves in such a direction 
so az to be contained within the angle defined by the angle of emer- 


gence vectors. 


It appears that there is always a zero position angle between 0 
and 90 degrees (the first quadrant with respect to the complex pole) 
such that the root locus from the complex pole closes on the complex 
zero. Thus a self-adaptive zero compensator could be located within 
this 90 degree sector and effectively "compensate" the poles. The 
residues from the roots of these closed root loci would be negligible 


compared with the dominant roots of the control system. 
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APPENDIX I 


Two digital computer programs were used. The first was "Program 
ANGMERG" to determine the angle of emergence for various pole and 
zero locations. Figures 7 and 8 were outputs of this program which 
verified analytical developments of Chapter II. 

The second digital computer program was ''Program Eqn" which is 
a program designed to plot root locus for any given characteristic 
equation (up to order 30). This program was originally programmed by 
LT Joseph D. Fennick, Jr., USN. A few changes in the original sub- 
routine were made to facilitate data input, reduce printout, and speed 
up graph plotting time. Figures 11, 12, and 13 are sample outputs of 
this program. Table 2 contains all of the input data (characteristic 


equations) for each of the 215 individual cases studied. 


63 


ee JOBN4GI4F sREESs EeGe 

PROGRAM ANGMERG 

THIS PROGRAM CALCULATES THE ANGLE OF EMERGENCE FOR VARIOUS 
POLE AND ZERO LOCATIONSe MULTIPLE RUNS CAN BE MADE BY 
INDICATING NUMRFR OF RUNS AND ADDING ADDITIONAL DATA CARDS 
REAL AN®) IMAGINERY PARTS OF POLES AND ZEROS ARE ‘READ IN 

WITH STIPULATION THAT REP(1) AND YMP(1) BE COMPLEX POLE ABOUT 

WHICH THE ANGLE OF EMERGENCE IS CALCULATEDe COMPENSATION 
ZEROS ARE NOT READ IN» BUT ARE GENERATED IN THE PROGRAM.) 

EET (REP( 2) AND YMP42) (BE CONJUGATE OF (REPU AND YMG) 

PART OF THE PRINT OUT AND GRAPHS CAN BE SUPPRESSED BY 
CHANGING VALUES OF W AND V RESPECTIVELYe FOR A CURVE OF 
RADIUS VS ANGLE OF EMERGENCE AT A PARTICULAR ANGLE OF 
COMPENSATION ZERO FROM COMPLEX POLEs ANGLE AAG IS READ IN 

IN RADIANSe TWO GRAPHS OF ANGLE OF EMERGENCE VS ANGLE 
POSITION OF COMPENSATING ZEROS AT A GIVEN RADIUS IS AVAILABLEe 
RR1 AND RR2 ARE THE RADII FOR THESE CURVESe ADDITIONAL 
SYMBOLS ARE== CE TROIDs THE INTERSECTION OF @ASYMPTOTESs«, 
THETA -- ANGLE OF ASYMPTOTES FROM HORIZONTAL WHEN COM PENSATION 
ZEROS ARE ADDEDss ZETA -— ANGLE OF ASYMPTOTES FROM HORIZONTAL 

WITHOUT COMPENSATING ZEROS»ss ANGZ -— SUMMATION OF ANGLES TO 
ALL ZEROS FROM COMPLEX POLE s, ANGP = SUMMATION OF ANGLES TO 

ALL VPOLES (FROM, COMPLEX POLE. 

DATA READ IN AS FOLLOWS cccece 


NOS IZ NUMBER OF RUNS 

M»N Zi2 NOs OF POLFS Ms NOw OF NONe COMP. ZEROS N 
ert) BF 1063 REAL PART OF POLES 
YMP(T) SF 10.3 IMAGe PART OF POLES 
REZ Cl) SELO63 REAL PART OF “<ZEROS 
YMZ(1) 8F 103 IMAG. PART OF ZEROS 

RR1»sRR2s AAGs Ws V 5F10.e3> RRI AND RR2 DESIRED RADIUS EO 


GRAPHs AAG IS ANGLE IN RADIANS LOCATING POSITION OF 
COMPENSATING ZERO FOR RADIUS VS ANGLE OF EMERGENCE GRAPH 
W SUPPRESSES PRINTOUT AS FOLLWS (W = 0¢0 PRINTS ONLY EIGHT 
VALUESs W = 1460 PRINTS OUT SIXTEEN VALUES FOR EACH RADIUS) 
V SUPPRESSES GRAPHS AS FOLLOWS ( V = 040 ALL GRAPHS 
SUPPRESSED» V = 1e0 A TOTAL OF THREE GRAPHS PER RUN 
DIMENSION REP(8) »YMP(8)sREZ(8) sYMZ(8) sR(900) 98(900) »REZ1(900) 
1ZANG(900) »ANG(900) sITITLE(12)9YMZ1(900) »ZA(8) 9PA(8) sGNA(900) » 
2 N(9NN)s AZ(9NN)» CEN(9N0)s ACF(900) 
READ 639 NOS 
63 FORMAT(I2) 
DO 492 K= 1sNOS 
READ 29MoeN 
2° FORMAT. (212) 
READ 45 (REP(I)9I = 19M) 
RFAD Gs (YMP(I)soI= 19M) 
IF(N=-1) 499599 
9 RFAD 459 (REZ(1)eI=1l9N) 
‘-READ G» (YMZ(1)9 I= 129N) 
4 FORMAT (8F19+3) 
READ 44% RR1s RR2» AAGosWoV 
44 FORMAT (5F1003) ° 


G PRINT DATA READ IN 
PRINT 629K 


CVV Y OY CVO EV CV ON CCV OY CVC OY OVO) CVO) 6) OV. OO OO ree re tion 


oly 


62 FORMAT (/4X sT7THRUN NOosI2) 
PRINT 69M N 
6 FORMAT (//9X52HM=513920X%s2HN=913) 
PRINT 8 
8 FORMAT (/5Xs7H REP(1)s9Xs7H YMP(1)59X97H REZ(1)99Xs7H YMZ(1)) 
IF(N=-1) 48549549 
Cee tN) 99s (REPCT Is YMP(1 is. T= 14M) 
BOTT 99 
49 PPINT 105 (REP(I)sYMP(I)>REZ( 1) sYMZ(1) 5] =19N) 
Ea eee] 
PRINT 99s (REP(I)s YMP(I)s IT=LLsM) 
OO UPORMAT (2F) S53) 
LO -PORMAT C4F15.3) 
PRINT 455 RR1sRR2sAAGoWoV 
45 FORMAT (/74Xs5HRR] =9F 10039 SXsSHRR2 =9F100359 SXs5HAAG =5F1003s 
15Xs 3HW =59F106355Xs3HV =5F 1003) 
XX = M=-N 
Kate AK = 2 
ASYMPTOTE VFCTOR WITH COMPENSATING ZFROS,5 THETA 
Pe. CX) 220252025. 203 
202 POINT 205 
205 FORMAT(/2X%s15HTHFTA= INFINITY) 
G60 TO 206 
203 THFTA = 180./X 
PRINT -2045s THETA 
204 FORMAT (/2Xs6HTHFTA=s F10-3) 
ASYPTOTE- “VECTOR WITHOUT COMPENSATING ZEROS» ZETA 
ZOG Ir “(M=aN) 210s 2105 211 
ZO PRINT’ 272 
212 FORMAT (/2X514HZFTA= INFINITY) 
GO TO 240 
21 il ZETA = 1800 /Xk% 
PRINT 214, ZETA 
214 FORMAT (/2X ,5HZETA=5F10¢3) 
ASYMPTOTE CFNTROIN»« CEN = (SUM =POLFS — SUM ZFROS)/ M-N 
240 SUM=0.0 
DO =201l=ls8 
201 SUM=SUM + REP(IT)-REZ(1) 
; SMM=-SUM 
ANGLES DUE TO NON-COMP. ZEROS, ANGZ 
ANGZ=0.0 
LEON ee SOs 250s 251) 
2G OUl2 1. = LEN 
PE CYMZ C1 =YMP UL) isi 2673 
71 °“TFCREZ(TI)=REPULIV 24s 76577 
74-201 VEATANF (LYMP(1)=YMZ01) ) 7 CREP] J=REZCI)))43624159 
GOTO iz 
te LPCREZ(T)=REP (1 )) 75380 680 
15> ZACT 123614159 
GO-TO 2 
76 ZAC(T 23 e/2e%32e14159 
GO 1G AZ 
74 2A) =20% 30 14 159—ATANF CUCYMPULj=YMZ 11 ))7 (REZ(1)=REPI1))) 
GOTO 12 
8O ZA(1)=0-0 
ea FA ia 
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13 TECREZ( by=REPIT)) 81582583 
Bl ZACIJ=ATANF(CYMZ( 1I=YMP (C1) )7IRE ZT) REP I) 


GO~70-"12 
82 ZA(1)=3414159/2. 
GO TO 12 
83 ZA(T) = ATANF((YMZ(1)-YMP(1))/(REZ(1) - REP(1))) 


12 -ANGZ=ANGZ+ZA(T) 
ANGLES DUE TO POLESs ANGP 
250 ANPG = 0-20 . 
IF (M=2) 2545 2545 253 
2534 PO 14 |] = 36M 
PPGYMe CTI =YMPUl) OTs 2393" 
91 IFCREP(I)-REP(1))94596397 
94 PA(T)=ATANF ((YMP(1L)-YMP(1))/(REP(1)-REP(1I)))4+3214159 
GO TO 14 
O72 FE CREP TL) HREP Ch) 9 551004100 
95 PA(1)=3214159 
GO TO 14 
96 PA(1)=30/2e%*3214159 
GO-TO 14 
97 PA(T)=20%*3e¢14159-ATANF ((YMP(1)-YMP(1))/(REP(I1)-REP(1))) 
GO TO 14 
100 PA(IT)=0.20 
GO TO 14 
93 IF(REP( I )—REP(1))1015102+103 
101 PA(I) = 3414159 = ATANF((YMP(1T)-YMP(1))/(REP(1) -— REP(I))) 
GO TO 14 
102 PA(1)=34e14159/2. 
GO TO 14 
103 PA(T) = ATANFCCYMP(T) =YMP(1))7(REP(1)—REP(1))) 
14 ANPG = ANPG + PA(T) 


254 ANGP = 3614159/2. 
ANGLE OF EMERG. IF COMPLEX ZEROS NOT PRESENT» ANGLE 
252 ANGLE=1806/3214159*(-3014159+ANGP—ANGZ ) 
GENERATE COMPENSATION ZEROS 
ALTe STATEs 15-19 FOR VARIABLE RADIUS 
DO 301 J=1s%20 
Cee 
R(J)=CA/10. 
15 DO 26 I=15630 


16 Al=I-1 

17 B(I)=AI/100.6 

Po oREZUC1 T= (REP CL). = COSe (Bt 1) iment) 
io YM ZC = YMP( 1) SINE CGT) eR) 


ANGLES DUE TO COMPENSATION ZEROS» ZANG(1) 

53 IF(B(1)-3614159) 5454955 

54 D(1)=1800 - B(1)*1806/3614159 
GO TO 56 | 

55 'D(1)=5406-B( 1) *1800/3014159 

56 IF(REZ1(1)-REP(1))57958959 

57 “ZANG(1)=1806/3614159*(3014159-ATANF((-YMZ1(1)-YMP(1))/ 
L(REP(1)-REZ1(1))))4+D(1) | 


GOSTO 20 
BG ZANGEE)=900e+BtT) 


GOTO .20 
59 ZANG(1)=180./3414159*(ATANF( (-YMZ1(1)-YMP(1))/(REZ1(1)-REP(1)))) 
Tot Dl) 
ANGLE OF EMERGENCF WITH COMPENSATION» ANG(T) 
29 GNA(T)=ANGLF~ZANG(1) 
ee Oy 22952295230 
230 CFN(I) = (SMM +2. * REZ1(1) )/¥ 
C2 omen) =S780. 7 3.14159 * (6s28 — B11)) 
2 UE UGNACIT)) 22323523 
22 ANG(I) = S=GNAI1) 
IF (ANG(I)) 60033335333 
609 ANG(I) = 360. + ANG(T) 
GO TO 333 
2ARANGIIY = 960.) — GNACT) 
Bao F CBC 1) RAG) 26456526 
46 ACF(J)= ANG(T) 
26 CONTINUE 
IF(R(J)-RR1) 3215279320 
320 JTF(R(J)-RR2)321527532]1 
27 LA=4HANG 
TE UVY 31753294531 
29 GO. 10 321 
2) DOT) Y=lsi12 
1 ITITLE(I)=8H 
ITITLF(1)=8H REFS 
ITITLE(2)=8HE.G. 
ITI TLE (3)=8HANGLF 
ITITLE(4)=8H OF 
ITITLF(5)=8H EMERG 
ITITLE(6)= 8H VS ZERO 
ITITLE (7) = 8H POSIT 
IF(R(J)-RR1) 33053315330 
331 ITITLE(8) = BH R=.5 
GO T6332 
330 ITITLE(8) = 8H R=1.0 
332 ITITLE(12)= 8H THESIS 
CAEL DRAW (63A0s AZ SANG sOs0sLAsTTITULE.,O0s0s2Z23 15 
Peis Oe FsSs5 Ce LAST) 
321 PRINT 405, R(J) 
405 FORMAT (/2Xs8HRADIUS =s5 Fl%e3) 
220 PRINT 28 
28 FORMAT (8X%614H ANGLE TO ZERO+s5xXs 19H ANGLE OF EMERGENCE, 5X 
1 9H CENTROID) 
IF(W) 60961460 
60 PRINT 305 AZ(1)sANG(1)9 CEN(1) 
PRINT 305 AZ(39) » ANG(39)59 CEN(39) 
PRINT 30+ AZ(78) » ANG (78) » CEN(78) 
PRINT 30s AZ(117)% ANG(117)s CEN(117) 
PRINT o204 AZ(157). ANG €157)4-°CEN (157) 
PRINT 305 AZ(196).s ANG(196)s CFEN(196) 
PRINT 3% AZ2(235)9 ANG(235)5 CEN(235) 
PRINT 309 AZ(276)s ANG(276)9 CFN(276) 
PRINT 309 AZ(314)5 ANG(314)5 CEN(314) 
PRINT 309 AZ(353)59 ANG(353)9 CEN(353) 
PRINT 305 AZ (392)% ANG(392)5 CEN(392) 
PAINT 426s AZ(423118 ANGI423)5 CEN( 431) 
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PRINT 30s AZ(471)9 ANG(471)>5 CEN(471) 
PRINE 30s JAZ (5101}s- ANG(510\s CENCS 107) 
PRINT 30s AZ(549)5 ANG (549)5 CFN(549) 
PRINT 30s AZ(588)5 ANG(588)>5 CEN( 588) 
GO TO 30 
61 PRINT 3205 AZ(1)sANG(1)9 CEN(1) 
PRINT 305 AZ(78) » ANG (78) »s CEN(78) 
PRINT 305 AZ(157)s “ANG 1167), CEN. (857 9" 
PRINT 30s AZ(235) 5s ANG (235) 5-CFN(235) 
PRINT 30s AZ(314)s ANG(314)s CEN( 314) 
PRINT 309 AZ (392)s ANG(392)s CEN( 392) 
PRINT 30s AZ2(471)s ANG(471)5 CEN(471) 
PRINT 309 AZ(549)s ANG (549)5 CEN( 549) 
30 FORMAT (3F20-3) 
30:1 (CONT INVUE 
PRINT 4015 ANGZs ANGP 
401 FORMAT(/6HANGZ =3F10¢395Xs6HAMGP =5F1003) 
PRINT 329 ANGLE 
32 FORMAT (/26HANGLE OF EMERGe W/O COMP =9F10¢3) 
IF V = 060s NO GRAPHe IF V = 1¢0 A GRAPH OF R VS ANG OF EMERG 
IF(V) 66365366 ‘ 
65° GO TO 402 
66 LA=4HANG 
LA =4HANG 
DO 47 I = 1l9l2 
47 ITITLE (1)=8H 
LTIVUECIL)S=8H -REES 
ITITEE( 2) =8HE «Gs 
ITITLE(3)= 8B8HRADIUS 
ITITLE(41= 8H- VS 
ITITLECS) = BH ANGLE 
ITTTVECEC6G) =. SHOF EMERG 
CALL DRAW (203sRo ACEsO0s0sLAsITITLEs090s25l9 
10s09799s0sLAST) 
402 CONTINUE 


67 END 
END 
il 
ly 
1260 1.0 1269 
~629 620 
eZ 204 —-3e7T ilo) 
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ee JORBRN4GI4GF, REFS Fe. Ge 


240 


24) 


206 
68 
70 
one) 

209 

203 

204 

324 


28 


205 
207 
22 


23 


20) 
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432 


PROGRAM EQN 

READ 2405 NUMB 

FORMAT (12) 

DO 241 I = 1» NUMB 

CALL SPLANE 

FND 

SUBPOUTINE SPLANF 

CI MENS TONNE (129% (129 )'s 11 (1014 ROOTR( 128)5RO0OTI (128). ITITUETI2(. 


POC 1271 sh C129) sPOOTI( 126) s5ROOTM( 128 14.AP(129)6AZ1129) 


COMMON ReVAR»sNO,ROOTRSZROOTI »sAP 4X 
MOND=1 

LAR =4H 

FORMAT(6E15.5) 

FORMAT (7+ 1l14HIMAGINARY PART,//) 
FORMAT (/7/s9HRFAL PART >s//) 
FORMAT (6AB) 

FORMAT(13) 

FORMAT(F1N.5) 

PEMD 7 OOst TITLE tT ye leleS) 

PRINT 324 

FORMAT ( JH]1s1l1HGRAPH TITLFs/) 
Pe Com in Lreetb ys l=156) 
FORMAT (/»36HORDER OF THE CHARACTERISTIC EQUATIONs/) 
PRINT 28 

READ 203sNO 

PRINT 2035NO 

N=NO+41 

FORMAT (8E19.5) 

FORMAT CBE A? oo.) 

FORMAT (is 53RCONSTANT PART OF THE COEFFICIENTS IN DESCENDING GROE 


IR»s/) 


PRINT 22 

REO e727 ) Se ACK ba Kad en) 

PRINT 207s5(A(K) 9K=1l9N) 

FORMAT (7s S3RVARTAGLE PART OF THE COEFFICIENTS IN DESCENDING: ORDE 


et Ae 


PRINT <2 3 

READ 205+(B(K)+K=1l +N) 
PRINT 207,(B(K) »K=l oN) 
VAR = 2001 

ND = 8 

XSCALF = 2.0 
YSCALE=SXSCALE 

FORMAT (/77/77s21HOTHE SYSTEM POLES: ARE. 74 
PRINT 201 

M=N 

DO 67 K=lo9N 

AP(K)=A(M) 

M=M—1] 


.- DO 432 K=1+530 


ROOTI(K)=06 


ROOTR(K)=O06 
NX=NO 

DO 15 K=143 1 
R(K)SA.A 


CALL POLYRT (APsXsNXsROOTR  sROOTI 51 e¢E—5) 
PRINT 69 ; 
PRINT 68s(ROOTR(K) sK=19NO) 
PRINT 70 
PRINT 689(ROOTI(K) »sK=15NO) 
CALL DRAW(NOsROOTR »sROOTI »MOD91 »LABsITITLEsXSCALE sYSCALE » 
lls7s2s29899s0sLAST) 
MOD=2 
202 FORMAT (/////7/521H THE SYSTEM ZEROS ARE>/) 
K=] 
Ss LPCBtk) > le2s) 
2 Keak +d 
GO TO 3 
1 NORD=N-K 
IF (NORD=-1) 69495 
4 ZFRO=-B(K+1)/B(K) 
7° FORMAT 7s 16HTHE SYSTEM ZERO=sF 1525, ) 
PRINT 7sZERO 
ROOTM(1)=ZERO 
ROOTM(2)=16e*XSCALE 
ROOTJ(1)=0-6 
~ROOTI(2)=06 
NORN=2 
GOT Os 11 
6 PRINT 9 
9 FORMAT (/+s25HALL ZEROS ARE AT INFINITY} 
GOr.10' 8 
5 NN=NORD+4+1 
DO 10 L=lsNN 


R(E)=BCK) 
10 K=K+1 
PRING 5202 
M=NN 
DO 46 K=1,5NN 
AZ(K)=R(M) 
46 M=M—] 


DO 9433 >K=h330 
ROOTM(K)=06 
433 ROOTJ(K)=O6 

NORX=NORD 
DO L6 Kel s3 1 

16 X(K)=0.0 
CALL POLYRT(AZsXsNORX sROOTM,ROOTJ510e¢E-5) 
PRINT 69 
PRINT 689{(ROOTM(K) »K=1sNORD) 
PRINT 70 
PRINT 685(ROOTJ(K) »K=1 5NORD) 

11 CALL DRAW(NORDsROOTM, ROOT OD Ener Ee ar ar YSCALE»s 
Lis7s2929839509LAST) 

8 CONTINUE 

GO T0(3193259335345 3593659375 38539540) »ND 

212G=160076 
GO 70 "41 

32 G=12016 
GO Topo) 

34a REis6pAB 
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GOr TO 4) 
Ba G=140312 
GO2TO 4) 
35 G=1.9394 
GO TO G4) 
36 G=1.9483 
GO) 17a.) 
GO TO 24) 
38 G=1.0633 
GO. TO 04s 
39 -G=1.071 
(ase ac aves 
40 G=1.97B 
41 CONTINUF ; 
DOr Ol ga 1435 
GO TOO" K= 15510 
DO. 300 L=l«n 
300 R(L)=A(L)+B(L)*VAR 
CAG ROOTS 
DO47 1 IIe lis. NO 
IF (ABSF (ROOTI(JJ))-52E-04) 61561562 
Gl vROOTI( Jo y=0. 
GOTO: 7} 
62 ROOTJ( JJ) =ROOTI (JUS 
71 CONTINUE 
100 VAR=VAR#G 
LE(J=-30) 101598498 
98 MON=3 
10) CALL DRAW(NOsROOTR,ROOTJ»sMODs2,LABs ITITLEsXSCALEsYSCALE y 
Lle7%s29298+s9s0sLAST) 
GO T0 206 
END 
SUBROUTINE ROOTX 
DIMENSION €(31)2D(29) »R(129) sROOTR(128) sROOTI(128) 5EE(31)sFE(31) 
PeAPCI29 lak 129) 
COMMON RoeVARsNOsROOTRsROOTI sAP 4X 


nee 
DO 1492 MA=1 43) 
ee CMaAy =(). 


1492 FF(MA)=0. 
20 BETAN=ROOTI(M)+EF(M) 
ALFAN=ROOTR(M)+FE(M) 
DO 7 =a loo 


S=2e*ALFAN 
T=-—(ALFAN#®*2+BE TAN#*# 2) 
CAI VERA 19 
C(2)=R(2)+S5*R(1) 
NC =NO+1] 
DO 2 L=25NC 
2 OG) PSR SEC CCK] De TEC VLD) 


AN= C(NO+1)-ALFAN*C(NO) 
BN= BETAN¥C(NO) 
IF (NO-3) 219 175 18 
17 CN 3e*R(1)*(ALFAN¥*2-BETAN¥*2) + 20*R(2)*ALFAN + R(3) 
BN Be FR(LIHALFANHBETAN # 20*R(2)*BETAN 


mou 


71 


21 


Is 


Ww 


19 


25 
bao) 


Orit 


434 


GO .TO E19 

CN = 2e*R(1)*ALFAN + R(2) 

DN = 2e*R(1)*BETAN 

GO TO 19 

D{ 1) = ct )) 

D(2)=C(2)+S*D(1) 

NU=NO=—1 

DO 3 N=39NU 

DIN) =C(N)+S*D(N-1)+T#D(N-2) 

CN= C(NO)—-22*D(NO-2) *BETAN*#2 

DN= 2e¢*BETAN* (D(NO-1)-ALFAN*D(NO-2) ) 

ALFA=ALFAN—(AN#*CN+BN*#DN) /(CN##24+DN*¥*2 ) 

BETA =BETAN + (AN*®DN=BN#CN)Z(CN*¥*#24+DN* #2 ) 

AZETA=e01+ABSFI(BETAN) 

ZETA=SORTF(201=T } 

IF(ABSF((ALFA-ALFAN) /ZETA )—e1)55955952 
“IF(ABSF((BETA-BETAN) /AZETA)=le) 56956952 

EE(M) =(ALFA-ALFAN)*10¢ + BETA -— ROOTI(M) 

FE(M)=(BETA — BETAN)*¥10¢ +ALFA = ROOTR(M) 

IF (A8SF((ALFA=ALFAN) /(ALFAN41.))—-5eE-4) 49495 

IF(ABSF((BETA-BETAN) /(BETAN+1.))~-5eE-4)69695 

ALFAN=ALFA 

BETAN=BETA 

CONTINUE 

MR=NC 

DO 6777 K=loeNC 

AP(K)=R(MR) 

MR=MR-1 

DO 434 K=1530 

ROOTI(K)=O6« 

ROOTR(K)=06 

NX=NO 

DO 83 K=1l»s31 

X(K)=0.0 

CALL POLYRT (AP»XsNXsROOTRsROOTIs1leE-5) 

GO T0 16 

ROOTR(M)=ALFA 

ROOTI(M)=BETA 

IF(M=$NO) 155916916 

M=M+1 

GO TO 20 

RETURN 

END 

END 


EeGe TRESTS 


19.0 Lil 2136 1488. 1088. 
1.0 320 18.225 


6 
10 


24e 22l1e 978. 28880 65286 
1.0 10 


72 


16625 


a 
- + 
+ 
if i : 
4 + 
+ 
+- 
+ ; ‘ 
5 re3a Be ‘ 
+ i 
=r “te 
+ ! 
“e 
= 
: a4) 
& 
Cs ptt mt 1 meer er —— cert errr te A 
~ OG eo lee on “We “39 
+- 
wt | 
= 
ne 
rome 
+ T 
+ a 
a n 2“ 
+ Ts 2S 
+ 
+ 2 
+ + 
= re 
+ F-3b + 8 
— 
* ~ 
- 
tT 
+ S&S 
T © 
eng ee 6S =O ~552 
> 
“- 


Fig. 11 Root Locus Data (Computer) 


73 


+ 

+ 
+ 
+ 
+ 
+ 
T 
Mi F=3c 
a 


-t 


+44 


(0074 
v. ; 


“ 


- 

+ 

+ 

aie 

+ 4 = 
< eee 
u 

2 

+ 

+ 

+ 


Fig. 12 Root Locus Data (Computer) 


74 


C52 


4 +h 
; 
Pa 


G2 








<i 





TABLE 2 


CHARACTERISTIC EQUATIONS 





i] 8 i] ‘ ‘ 
i] i] t L] § 
Q omonworvso i WWMWWMW ! momo nwowMmrsd i! 0 0 0 wy t oa oO ™~ 
WNAN OO St ONNNN WONMN ANON St WNNOHAON On Ono 
t * 28 e@© © 2© § © e+ @& @ euers 7. 2©«© © @© «© #@® # @ 9 eo e« e@ © e@ © @ Q oo @ e@ o e@ 
eat OMWO oon wom Aon Tr KH OMOM OM Oo OrRAAWY WY ONO + © 
4 NN Fee testa Ness tC Nes ANN Fees (OO Ftxt Xi ketiAtesst IP OANNN =| 
! —. os perwee A mem pri enrerrmerenmecnence arrencere Avene mee renee ts RAN Toran SN RAS OR weonroent 2 er res (oe re etadien ahi de) bce iatean lott “ET. AIS Nae FE Hae! 
i] ‘ t] t § 
N wo 3s 4 § ( N co N 4 
co 0G ° foe) CON foe) oO N foe) oO N onmvo wn oO Ke) foe) 
Om ANN 1 e °* © § © ' ©€ ° e ° * t @©M © © «© @ °° § € tt 
© aeGgANnNN ON tAIN ON DOF ANN MON AN I OFA ANNNN ONOMONON 
=i ‘ ie ‘ert ‘(ee ire 
ce RCI TP Et NAT NE RETRO SIE SES AE LD TOE TE DAE AAS IE EL IE OL ES ETOH A PLO RE LE ELODIE DTT T Srey i oe Ta REPT PIE owns air pie eee Cas at or 
§ t] t] L} 
4 ‘ t a 
Ln oe) § Oo 'o 1 © 
rw ae) it eos 
| on aeeetenenni eh eieeiaenieal ied ~ = pera trect ane i are hes eee tae Tear sa a Q Te f] 
‘ 4 4 t- 
(ie tN a 90} a 8) 
ise) + ce ion! 
(ere ior te tee 
I et CNG Ae CITE AEE Ay RE IC Tt DES CE EWE EV ST ee haan PAGANS UT a + see A mia 4 orrebB a | edna oa = anh 
i] t L) t 
i) ‘© ' © 1 © 
pS Gee ag ee ee eee Nn ene arena Sapna: — Le eva nee ae, 
‘ t t ‘ 
‘ i] 4 t 
ee i a : eae _ 
~~ meneeas creer raver eres rents emer er rere fre eer cerns a a scar inate 
‘ 4 ‘ ‘ 
‘ ‘ 4 J : ” _ 
<i : nectar : j 
OnonWnoWNW ! oooo°o !} Owmnownono vr ‘ OoOwnown oe 
Oot DMOMON OoOneonm OTFDMNONKH I ot OM © 
Che Sa ORO A Lt 6 Oo Sr eS * '*oOoOne 














——— Te Wrage eae cael AER LB ALE A EE LES 





i] 

(] 
GO 0T DH OD e020 T SCOUT VH 00.4 GOUT DH OO GC O0T89 VY 
Sn Fe rearrested a BS aL A eA A I a: a ea Sa Cie: alae an Atel ae aula ell ea eb aria ary) 
fdeaq edaddeded eetdd  adddtddee tee detdd aadaas 





76 
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